Efficient Grover search with Rydberg blockade 
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We present efficient methods to implement the quantum computing Grover search algorithm using 
the Rydberg blockade interaction. We show that simple 7r-pulse excitation sequences between ground 
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different properties of the atomic interactions are identified and the error scaling of the protocols 
with system size is found to be promising for experimental investigation. 
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I. INTRODUCTION 

The last few years have witnessed impressive progress 
in attempts to perform quantum computing with trapped 
neutral atoms via the strong, long range Rydberg block- 
ade interactions [IHS], and research has begun to scale 
the physical systems to quantum registers encoded in 
several trapped atoms and to perform non-trivial algo- 
rithms. The Rydberg blockade gate was proposed in [6] 
as a robust and fast quantum gate, relying on the pres- 
ence of a single atom in a high lying Rydberg state shift- 
ing the Rydberg excited state energy of any other nearby 
atom and thus preventing the resonant excitation of that 
atom. As the first (control) atom can be selectively ex- 
cited from any of the ground states representing qubit 
states |0) and |1), the evolution of the second (target) 
atom is effectively controlled by the quantum state of 
the first one. 

With the capability to perform suitable one- and two- 
bit gates, one can decompose any quantum computation 
as a suitable sequence of such gates. The Rydberg block- 
ade interaction mechanism has a special property that 
we will make use of in the present communication: A sin- 
gle Rydberg excited atom can simultaneously block the 
excitation of a whole ensemble of atoms in its vicinity. 
This collective blockade has been suggested as a means 
to use an ensemble of atoms to code a single qubit in [7] , 
to generate a wider class of entangled multi-atom states 
[U [9] , and to code multi-bit quantum registers in ensem- 
bles of multi- level atoms [lOj [11] . In this article we shall 
show that also in the more conventional encoding with 
one qubit per atom, the simultaneous blockade among 
many atoms offers interesting prospects for unique multi- 
bit quantum gates. In particular, we propose means to 
carry out the Grover search algorithm [12] on a k-h\t 
register, encoded in k individually trapped atoms. 

The outline of the paper is as follows. In Sec. II, we dis- 
cuss the theoretical Grover search algorithm, we describe 
the Rydberg blockade gate mechanism, and we show for- 
mally how the two key operations in the Grover algorithm 
can both be implemented efficiently by sequences of tt- 



pulse excitations between the qubit levels and the Ryd- 
berg states. In Sec. Ill, we discuss possible alterations 
of the gate sequences suited to the case of atoms with 
different interaction properties. In Sec. IV, we estimate 
the expected fidelity of the gate operations proposed in 
the manuscript. Sec. V concludes the paper. 



II. GROVER ALGORITHM 

The Grover search algorithm is able to identify with 
high probability a single marked element xq out of N 
candidates with only queries of a database or a suit- 
able oracle device. The algorithm assumes that a register, 
holding a superposition of states ^ x c x \x) corresponding 
to the different possible values x of the unknown variable, 
is queried by being acted upon with a unitary operation 
that leaves all components unchanged, except \x$) which 
undergoes a change of sign, c XQ — >• — c Xo . The query pro- 
cess thus yields a conditional phase shift and does not 
change the weight of the desired state component within 
the superposition, and hence it does not change its prob- 
ability to reveal itself in a measurement on the register. 

The key ingredient in the Grover algorithm [12 is the 
subsequent step, which replaces all amplitudes c x by their 
values reflected in their mean value c = ^ x c x : 

° x ~^ ^ ^2° x ' ~ (° x ° x ^' w 

x' 

It is clear that if all amplitudes are initially identical with 
the value 1 / y/N, after the conditional change of sign and 
the reflection in the mean value, the magnitude of the 
marked element has changed. For a large value of AT, 
one gets c Xo ~ 3/y/N while all other amplitudes are still 
of the order 1/y/N after the first step of the algorithm, 
and repeated action of the steps will further increase the 
amplitudes of the marked element by of the order 1 / 
in each iteration. The \xo) amplitude approaches unity 
after ~ steps [T2] . 
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FIG. 1: Atomic configuration for quantum computing with 
Rydberg blockade gates, a) An ensemble of atoms individu- 
ally trapped with mutual distances smaller than the Rydberg 
interaction blockade radius, b) Level scheme of individual 
atoms, showing two stable states |0) and |1) used for qubit 
storage, and a Rydberg excited state \r) which may be ex- 
cited either directly or, more conveniently (dashed arrows), 
by a two-photon transition via the optically excited state \e). 
c) The blockade interaction is due to the the dipole-dipole in- 
teraction which shifts the Rydberg excitation energy of atom 
j if atom i is already excited. 



A. Query of the oracle with the Rydberg blockade 
interactions 

We now assume that the register is represented by 
k = log 2 N individually trapped atoms, see Figl.a), with 
two low-lying long lived states |0) and |1), and a Rydberg 
excited state |r), see Fig.l.b). The arguments x of the 
database search are given in a binary representation, and 
in particular the marked element xq has the binary rep- 
resentation 60,61, ...6fe_i, with bi = or 1, represented 
by atom i occupying state |0) and |1), respectively. To 
certify whether a candidate x equals xq one must there- 
fore check that every bit value in the representation of x 
coincides with the corresponding bi in x$. In an exper- 
imental implementation, we imagine that we hand over 
the physical system to the oracle - another component of 
the experiment, delegated to this role. The oracle must 
then perform the unitary operation, which puts a change 
of sign on the xq component and leaves all other state 
components unchanged. 

Here is how this can be done: the oracle component of 
the experiment is in possession of the values 60, 61, 1 
and can thus apply operations on the individual atoms, 
where the ground state amplitude in |1 — bi) = \0i) or 1 1^ ) 
is transferred by a resonant 7r-pulse to the Rydberg ex- 
cited state \ri) and back by a second 7r-pulse on the same 
transition. The accumulated 27r-pulse yields the initial 
state with a minus sign if the atom actually occupies the 
state |1 — bi) but no phase change if the atom occupies 
the opposite state, and this operation thus provides the 



conditional phase component of the Grover search on a 
one-bit register. 

We are interested in testing whether the input state 
simultaneously fulfills agreement with all bits in the com- 
parison and here the Rydberg blockade is useful. Fig.l.c) 
summarizes the basic idea: if the control atom i is excited 
to a Rydberg state, the excitation energy of the Rydberg 
state of the target atom j is raised by the strong dipole- 
dipole interaction between the atoms, and hence excita- 
tion by a laser which is resonant with the single target 
atom transition is blocked. A 2tt pulse on the target atom 
thus has two possible outcomes: no effect if the control 
atom occupies the Rydberg state and a change of sign if 
the control atom does not occupy the Rydberg state. 

Applying the 7r-pulse |1 — bo) — )> \tq) will result in ex- 
citation of the components of the register quantum su- 
perposition state where the th qubit does not match the 
marked element. One proceeds with a 7r-pulse |1 — b\) — >• 
1 7*1) on the I s * atom. Note that due to the blockade, this 
pulse has no effect on state vector components where the 
th atom is already excited, but components which did 
match the marked element on the th bit, and which were 
therefore not excited, will now become excited if they do 
not match the 1 st bit. This is illustrated in Fig. 2, where 
we use the symbol Ro^fjr) to denote a 7r-pulse excitation 
between states |0) and |r), and similarly for state |1). 
The figure shows how we attempt to successively excite 
each atom from the ideally unpopulated state |1 — bi) to 
the Rydberg state, and eventually all input state compo- 
nents have undergone a unitary transform into new states 
with precisely one Rydberg excited atom (the first atom 
in the register not fulfilling the comparison with the de- 
sired bit string) or no Rydberg excitation at all in case all 
atoms occupy the appropriate ground state \bi) defined 
by the marked element. Applying now a second set of 
7r-pulses in reverse order, starting with the k — 1 st atom 
and finishing with the th atom, we return all atoms to 
their initial state, and we equip all components that ac- 
tually underwent a 2tt transition to the Rydberg state 
and back with a change of sign. Note that if we had in- 
tertwined a resonant Raman process on an ancilla target 
atom between states |0) and |1) via the Rydberg state 
between the two sequences of 7r-pulses, the net effect of 
the operation would have been a /c-atom controlled NOT 
operation on that atom, controlled by the agreement of 
the fc-bit atomic register with the bit string 60? •••6fc— 1? 
originally proposed by Isenhower and Saffman [13]. 

Since the Rydberg blockade ensures that components 
of the quantum state with any number of register values 
differing from the b^s will be excited and de-excited at 
precisely one location (atom) in the register, while the 
\xo) component is left unchanged, the whole operation 
can be written 

l*> - < -fi), xoV (2) 

Apart from an irrelevant global change of sign, this ac- 
complishes the conditional sign step of the Grover algo- 
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FIG. 2: Grover conditional phase: Successive 7r-pulse excita- 
tion transfer of all atoms from one of their qubit states to the 
Rydberg state. The lasers couple the states |1 — 6$), comple- 
mentary to the bit values of the marked element xo, to the 
Rydberg states. If one atom is already excited, no further ex- 
citation occurs, and hence quantum register components with 
more atoms populating the states coupled to the laser field 
will only be excited once, while the state with no atoms cou- 
pled to the laser fields will not be excited at all. A second set 
of 7r pulses, applied to the atoms in reverse order, returns all 
population to the initial states, but a relative change of sign 
has been accumulated between the state \xq) and all other 
components. 



rithm. 



B. Inversion about the mean with Rydberg 
blockade interaction 

The inversion about the mean ([!]) is given by a sim- 
ple expression in the ^-representation. If we stay in the 
Hilbert space of dimension N — 2 fc , the transformation 
is written, |*) — >► E/g|*)> with the matrix 

U G = 2P- J, (3) 

where P is the N x N matrix with the value 1/N in all 
positions and / is the N x N identity matrix. Defin- 
ing the initial state of the search algorithm with uniform 
amplitudes, |^o) = J2 X ~^\ x )i we observe that P is the 
projection operator, 

P=|*oX*o|. (4) 

Note that we can write 

U G = P-{I-P)=P-Q, (5) 

where Q is the projection on the space orthogonal to 
|^o)- Writing the operation like this, we recognize that 
this is also, apart from a global phase factor, a conditional 
change of sign, conditioned on the state not being |^o)- 



The quantum algorithm is applied on an atomic quan- 
tum register, where different input values are encoded 
in binary representation in the individual ground state 
atomic manifolds. The uniform input state |^o) and the 
associated projection operator P is so far expressed in 
the basis of the variables x, but it is also straightforward 
to identify the expansion of the state in binary represen- 
tation, where |^o) is simply a product state of all qubits, 
prepared in the superposition state (|0) + \l))/y/2. To 
multiply a phase factor on the state component along 
|^o) is thus equivalent to performing this multiplication 
conditioned on all qubits being in the state ( |0> + 1 1>)/a/2. 
This is, indeed, equivalent to our task in the previous 
section, except that in that section, we checked the occu- 
pation of the states |1 — b{) complementary to the desired 
classical bit value while now, we must check for the 
occupation of a superposition state. 

Fig. 3 illustrates how in the second Grover step, exci- 
tation from qubit superposition states is driven by two 
fields in a A-configuration with well defined phases and 
amplitudes, while in the first Grover step, shown in Fig. 2, 
checking for the bit values b{ uses only a single field from 
the |1 — hi) state. It is interesting to note, that the n- 
pulses applied always occur (or do not occur) from the 
ground state manifold to the Rydberg state and back, 
and not between superpositions of these states. This 
implies that a variety of pulse sequence schemes and 
chirped adiabatic passage techniques can be applied to 
make these pulses robust against small disturbances of 
the system and variations in the field parameters, and 
this robustness also applies to qubit superposition states 
in the ground state manifold [13] . 

The projection of the full register state on |^o) van- 
ishes if just one atom i does not occupy the state 
(\0i) + |li))/V% and we can hence check the state bit 
by bit by exciting the Rydberg state with two laser fields 
in a A-configuration, for which (|0i) + \ li))/y/2 is a dark 
state, i.e., with two Rabi frequencies of same magnitude 
and opposite sign [14]. This laser configuration excites 
the bright state \Bi) = (|0;) - \U))/V2 into the Rydberg 
state and serves the same purpose as our excitation of the 
|1 — hi) register state in the first step of the algorithm. 
As above, we thus begin the sequence of operations by 
exciting the 0^ atom coherently from the bright state 
to the Rydberg excited state |r*o). Thereafter, we excite 
the 1 st , and subsequently all the following atoms on the 
same transition, noting that for any component in the full 
register superposition state the first atom in the bright 
state will be excited and will hereafter block excitation 
of any further atoms in the same state component. The 
net result of the sequence of excitation steps is that the 
|^o) component of the register is left unchanged, while 
all other components have precisely one Rydberg excited 
atom. Applying the same 7r-pulses, but now in opposite 
order to the atoms in the register, will return all Rydberg 
excited state population back to the ground state where 
it came from, but writing also a minus sign on the quan- 
tum state. So, after the pulse sequences, all population is 
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FIG. 3: Grover inversion- about- the- mean: With the A- 
transition laser excitation scheme shown, the (|0i) + \ li))/\/2 
dark state in each atom is uncoupled while a 7r-pulse exci- 
tation transfers the bright superposition qubit state \Bi) = 
(\0i) — to the Rydberg state. In a succession of such 

excitation pulses, if one atom is already excited, no further ex- 
citation occurs, and hence quantum register components with 
more atoms populating the bright states will only be excited 
once, while the state |^o) with no bright state atoms will not 
be excited at all. A second set of ty pulses, applied to the 
atoms in reverse order, returns all population to the initial 
states, and a relative change of sign has been accumulated 
between the state |^o) and all other components. 



restored into the ground state qubit space, but all com- 
ponents of the state vector except |^o) have undergone 
a change of sign. This is precisely the unitary operation 
needed for the Grover inversion-about-the-mean! 

It is worth pointing out, that by leaving the ground 
state qubit space and carrying out multiple operations 
on several atoms before finally returning to the qubit 
space again, the entire multi-bit gates described in this 
and the previous subsection are explicitly not formed as 
sequences of one- and two-qubit gates. The Rydberg 
blockade interaction between all qubits provides a short- 
cut to multi-bit operations, which may not be found in 
other quantum computing implementations, where other 
shortcuts may apply instead, e.g., via the centre-of-mass 
motional degree of freedom in quantum computing with 
a string of trapped ions [15] . Indeed, some of the ingre- 
dients in the above analysis of the Grover algorithm were 
already discussed and can be implemented with quite dif- 
ferent interaction steps in the ion trap computer [16] and 
in cavity QED based proposals where a single field mode 
couples to all qubits [T7] . 



III. ALTERNATIVE CONFIGURATIONS AND 
INTERACTION SCHEMES 

Favorable conditions for Rydberg blockade are associ- 
ated with the Forster resonance phenomenon occurring 



when there is a near degeneracy between the product 
state with two Rydberg excited atoms and another Ry- 
dberg product state [18-20 . None of the single electron 
states have mean dipole moments, but the non-vanishing 
dipole matrix elements between the different single atom 
states couple the two product states and split their en- 
ergy levels. The resonant Forster resonance mechanism 
gives rise to large energy shifts scaling as 1/R 3 with the 
interatomic distance, approaching 1/R 6 for very large 
distances, where the coupling gets smaller than the en- 
ergy difference between the two states. Being relatively 
isotropic for s states, the Forster mechanism is useful to 
provide coupling between any pair of atoms in a not too 
large system of particles. 



A. Simultaneous excitation of all atoms 

It is tempting to look for a protocol which uses the 
same excitation mechanisms of the individual atoms, but 
applies the coupling lasers at the same time to all atoms. 
In both steps of the Grover algorithm, there is a state 
component which is not excited at all, and this com- 
ponent would also be invariant under the simultaneous 
application of the laser fields to all atoms in the regis- 
ter. The state components which differ from the regis- 
ter state \xo) = \bobi...bk-i) (or from the uniform state 
|^o)) at exactly one qubit location, will be excited pre- 
cisely at that single location under the simultaneous in- 
teraction with the laser fields, and after the subsequent 
deexcitation, the register acquires the desired phase fac- 
tor for these components. If driven simultaneously, the 
state with two "erroneous" qubits, z, j, experiences a cou- 
pling towards the state with two Rydberg excited atoms. 
That state, however, is shifted by the Forster interac- 
tion, and what happens, instead, is that the ground 
state, say \gi,gj) couples to a symmetric combination 
(\ r ii9j) + \9ii r j))/ \/2 with only one Rydberg excited 
atom. The Rabi frequencies of the two coupling terms 
interfere and lead to an effective coupling which is y/2 
larger than the single atom excitation Rabi frequency 
PP. If three atoms occupy states, simultaneously coupled 
to the Rydberg states, the enhancement is \/3, etc, and 
since the superposition state of the register contains com- 
ponents with all these different numbers, it is not possi- 
ble to apply a simple constant laser pulse, which will act 
as a 7r pulse on all the transitions. One may imagine 
a pulse sequence, identified by optimal control theory, 
which by suitably varying the amplitudes and phases of 
the fields, yields an effective 7r-pulse on all or many of 
the atomic components with different Rabi frequencies. 
When designing such a pulse sequence one must ensure 
that it does not widen the resonance band width beyond 
the Rydberg blockade interaction shift and hence permits 
doubly occupied Rydberg states. Sequential addressing 
of the individual atoms may ultimately be both faster 
and more robust than such a collective stategy. 

A more realistic scheme including simultaneous exci- 
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tation of all register atoms can be engineered if excita- 
tion is possible to two different Rydberg states, \s) and 
|r), with the property that a pair of state \s) atoms is 
not near degeneracy with any other pair of states (no 
Forster resonance), while a pair occupying \r) and \s) 
does constitute a dipole-dipole resonant state and hence 
experiences a large blockade shift. Possible states in ru- 
bidium with this property are \s) = 1 40p 3 / 2 i m = 1/2) 
and |r) = |41s 1/2 ,ra= 1/2) [9]. 

Some calculations indicate that the dipole interaction 
shifts and hence the blockade mechanism are generally 
present also for more than two excited atoms [21 , but 
special cases have been identified where two atoms block, 
while among three atoms resonant excitation of more 
than a single atom may be allowed [22]. We will assume 
that one can avoid the latter phenomenon by appropri- 
ate choice of atomic systems, but we note that it may 
take a more detailed calculation to precisely assess the 
collective Rydberg interaction shifts within a collection 
of state \s) excited atoms and a single \r) atom. 

Now, in this situation one can simultaneously drive 
7r-pulse excitation pulses on all atoms from the relevant 
single qubit states to the Rydberg state |s), which thus 
gets a total occupancy of zero for the state component 
\xq) and |^o) m the first and second part of the Grover 
step, respectively, and a total occupancy of unity or more 
(!) for all other components. We assume that a single 
ancillary atom is present in the state |0) and has not yet 
been excited by laser fields. This atom is now driven 
towards the Rydberg state |r), and due to the Forster 
resonance interaction, this transition is blocked if there 
is already any number of atoms (> 1) present in state 
see Fig. 4. 

Completing a full 2tt pulse on the |0) — \r) transition of 
the ancillary atom thus yields a change of sign if the regis- 
ter initially occupies the targeted state \xq) (first Grover 
step) or |^o) (second Grover step), and no process and 
no change of phase occurs otherwise. Transferring finally 
all atoms in state \s) simultaneously back to their ini- 
tial ground state by inversion of the first 7r-pulses thus 
completes the operation (note that we invert the 7r-pulses 
by laser fields with the opposite sign of the first exciting 
pulses - in this way, we do not accumulate sign changes 
of the state vector). A single iteration of the conditional 
phase shift and the inversion about the mean is thus ac- 
complished in the time it takes to perform four single 
atom 7r-pulses on the ground state to \s) transitions and 
two single atom 27r-pulses on the ground state to \r) tran- 
sition. In a database with N elements, one has to repeat 
these operations ~ y/N times to complete the full Grover 
search algorithm. 



B. Grover search with a sub-register architecture 

To scale the system to large numbers of qubit s, it is dif- 
ficult to retain individual addressability while still find- 
ing room for all the atoms within the Rydberg interaction 




FIG. 4: Grover conditional phase: Simultaneous 7r-pulse ex- 
citation transfer of all atoms from a specific qubit state to the 
Rydberg state \s). The \s) state atoms do not interact with 
each other, and the number of Rydberg excited atoms is not 
limited to zero and unity. Finally one ancilla atom, prepared 
in state |0) is excited by a 27r-pulse towards the Rydberg state 
|r), which does interact with all |s)-state atoms (see the in- 
sert sketch of the atomic mutual interactions). A second set 
of 7r pulses, applied to all atoms in reverse order, returns all 
atoms to their initial states, but a relative change of sign has 
been accumulated between the state with no Rydberg exci- 
tations and all other components. Exciting the qubits from 
the states complementary to the bit values bi thus yields the 
conditional phase operation, while excitation from the bright 
qubit superposition states \Bi) — (|0i) — \li))/y/2 yields the 
inversion about the mean. 



range of each other or of a central ancilla particle. In this 
section we will show, that it is then possible to split the 
system in separate sub-registers, which each constitute a 
cluster around an ancilla atom, as sketched in Fig. 4. As 
in Fig. 4, we assume no Forster resonance and hence no 
blockade effect associated with atoms in Rydberg state 
|s), while any non-zero number of logical qubit atoms ex- 
cited to \s) will block the excitation of the ancilla qubit 
to the Rydberg state |r). At each query and inversion 
step, we start with all the ancilla atoms in state |0), we 
excite the sub-register atoms to the state \s) as in the 
previous section, and we transfer the ancilla atoms via 
\r) to the state |1). The relevant property of the sub- 
ensemble logical qubits is thus encoded in the population 
of this final ancilla state rather than in a phase shift. The 
logical qubit atoms are then transferred back to their re- 
spective ground states by inverting the 7r-pulse excitation 
processes applied to each atom, and the task is now to 
compare the ancilla atoms. 

Changing the sign of the quantum state with all an- 
cillas occupying state |1) while retaining the sign on all 
other states would be easy if all ancillas were within the 
Rydberg blockade radius of each other: it would merely 
consist in applying the query operation, described in Sec. 
II. A with the ancilla register marked element Xq = 11... 1 
. This query of the ancillas could be obtained by physi- 
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cally moving the atoms from their sub-register locations 
to a tighter interaction volume, in the same spirit as 
the ion shuttling architecture [23 . Let us present, how- 
ever, a scheme that will work with stationary ancilla 
atoms which experience only a finite Rydberg interaction 
strength AE rr , which is too weak to exercise a reliable 
excitation blockade. With this interaction, the internal 
state of one atom can control a transfer process in an- 
other atom in the following way: i) the control atom is 
excited and the target atom is excited from a (bright) su- 
perposition (| a) — \b))/y/2 of the initial and desired final 
states | a) and \b) to the Rydberg state |r), ii) a change 
of sign accumulates on the \rr) product state component 
and Hi) inverse pulses return the Rydberg state ampli- 
tude to the ground states. The resulting controlled sign 
change on the bright superposition state of the target 
atom is equivalent to a transfer between the initial and 
final state [14J, for example, 



\a) = -L((\a) + \b))/V2 + (\a)-\b))/V2) 
-L((\a) + \b))/V2-(\a)-\b))/V2) = \b). 



(6) 



While not each ancilla atom needs to interact with all 
other ancilla atoms, we assume that we can pick a se- 
quence of operations involving always pairs with suffi- 
cient Rydberg interaction to accumulate a change of sign, 
and hence the controlled state transfer, after a dwell time 
r = 7r/AE rr in the Rydberg state of a few microseconds. 

We propose to verify the ancilla quantum states in a 
binary search tree pattern, so that we group the ancillas 
in pairs, encode their joint state into the first atom of 
each pair, form new pairs of these ancilla atoms for which 
the joint state is again encoded in the first atom, and 
repeat until one final atom occupies state |1) if and only if 
all the ancillas were put in state |1) by their sub-register 
qubits. A 2tt pulse on the last ancilla atom provides 
the desired conditional phase, and reversing all ancilla 
operations through all branches of the search tree, we 
return all ancillas to the initial state |0), retaining the 
conditional phase factor only on the marked element in 
the query operation and on the corresponding component 
in the inversion about the mean operation. 

Two ancilla atoms may populate four logical states, 
1 00), 1 01), 1 10), and |11), for which we want the first atom 
to end up in state |0), |0), |0), and |1), respectively. Since 
we want to apply reversible and unitary operations, for 
three product states to end up with the same state of the 
first ancilla atom, the second ancilla must accomodate 
three orthogonal states. We thus assume ancilla atoms 
with a three level structure |0), |1), |2), and note that an 
operation which transfers 1 10) to 1 02) and retains all other 
product states will be adequate. That operation can be 
done in two steps, 1 10) — >> |12) — >• |02), which are both 
of the controlled transfer type (6]): i) the state |1) of the 
first atom controls the transfer between |0) and |2) of the 
second atom, ii) the state |2) of the second atom controls 
the transfer between |1) and |0) of the first atom. By 



going through all four possible input states, one verifies 
that the operations indeed perform the desired task. 

As we set out to achieve, the first ancilla atom in each 
pair now encodes the logical AND of the two atoms, and 
by forming the new pairs and repeating the operations a 
total of n s /2+n s /4+n s /8... = n s — 1 times, a single atom 
witnesses the state of all n s sub-ensembles. We invert 
all steps by the same operations in reverse order, and 
hence a single query or inversion about the mean on the 
full register, involves 2 blockade gates between the qubit 
atoms and the ancilla atom within each sub-register (2n s 
operations carried out sequentially or in parallel), and 
2{n s — 1) interaction gates between pairs of ancilla atoms 
(which can also be done in fewer simultaneous operations 
on different pairs). 



IV. SCALING OF ERRORS 

The search size that can be implemented using multi- 
bit gates depends on the error scaling with k. A com- 
plete analysis of this question, even in the ideal situation 
of perfect hardware, would require a simulation of the 
multibit gate that accounts for errors due to spontaneous 
emission from Rydberg levels, finite blockade efficiency, 
and the multiplicity of high n Rydberg levels. Such an 
analysis is cumbersome and has not yet been performed, 
even for the much simpler case of a two qubit CNOT 
gate. We will instead follow the approach used in [5j [19] 
and assume the gate error averaged over all possible input 
states is small, so that the error can be well approximated 
by separately adding errors due to Rydberg state sponta- 
neous emission and finite blockade strength. When such 
an analysis is performed for a controlled phase gate it is 
found that the error is minimized by choosing a Rabi fre- 
quency ^ op t = (77r) 1 / 3 B 2 ^ 3 /r 1 / 3 which gives a minimum 
gate error E = (3(77r) 2 / 3 /8)(Br)~ 2 / 3 . Here B is the two- 
atom blockade interaction strength and r is the Rydberg 
spontaneous lifetime. 

Applying the same type of analysis to the Grover 
phase inversion step leads to equivalent scaling ^ opt ~ 
g2/3^ r i/3 anc j £ ~ c(/c)(Br) -2 / 3 with a numerical pref- 
actor c(k) which is approximately linear in the number 
of bits k. A detailed accounting of the analysis leading 
to this result will be presented elsewhere [24 j. In the 
present context it is interesting to note that there are 
qualitative differences in the error scaling of the multibit 
gates as compared to that of the two-bit CNOT gate, de- 
spite the mathematical similarity of the result. For the 
two-bit blockade gate intrinsic errors < 10 -4 are theoret- 
ically possible by exciting Rydberg levels with n > 100 
[25] at small atomic separations of R < 5 fim. Low errors 
can be obtained at larger separations by going to higher 
n. At large R we have a van der Waals interaction and 



B - n^/R 6 ,r - n 3 so E 



~ n 



-28/3^>4 



In contrast with this situation, implementation of a 
multibit interaction leads to some additional constraints. 
Let us suppose that k atoms are arranged on a two- 
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dimensional square lattice with y/k atoms on a side and 
period d. Ideally we want to have a strong inte raction fo r 
atoms with the maximal spacing of i? max = y/2(k — l)d. 
We may choose n large enough to ensure an adequate 
B(i^max)- However, doing so will result in B(d) being 
a factor (2(fc — l)) 3 larger for van der Waals scaling. 
If the resulting B(d) exceeds half of the spacing be- 
tween Rydberg levels of neighboring n then the block- 
ade can actually be reduced for proximal atoms, which 
will lead to large errors. To avoid this from happening 
we must respect B(d) < [U(n) — U(n — l)]/2 ~ n~ 3 . 
Thus the error scaling for a multibit gate is E ~ (n -3 x 
n 3 ) -2 / 3 ~ constant. The implication is that the prin- 
cipal quantum number n should be chosen large enough 
to ensure that the entire array is in the resonant dipole- 
dipole limit (B ~ l/^ 3 ) to minimize interaction strength 
variations across the array, but once we have met that 
constraint there is no advantage to going to higher n. 

We have analyzed the fidelity of a Grover phase in- 
version step, including lattice averaging of the different 
powers of the interaction strength which arise at inter- 
mediate stages in the error analysis. It turns out that it 
is not possible to limit the maximal interaction at spac- 
ing d and simultaneously have the entire lattice in the 
B ~ 1/R 3 regime for k > 9. We have therefore averaged 
the actual computed values of B over the lattice. The 
errors for different architectures are given in Table [TJ For 
Cs \ns) states with n = 75 we find that for search of a 512 
element database (k = 9) the average error is E ~ 0.002. 
The error for inversion about the mean will be essentially 
the same, so one Grover iteration has an expected error 
of E ~ 0.004. The error for 65536 elements (k = 16) 
is E ~ 0.015. It should be emphasized that these are 
average errors and the actual error for some input states 
may be much larger. 

The error scaling is different for the simultaneous ex- 



citation approach presented in Sec. |III A| In this case we 
require a large asymmetry between the \s) — \s) (weak) 
and | s) — \r) (strong) interactions. We choose states and 
lattice spacings such that \s) — \s) is in the van der Waals 
regime and scales as n 11 while \s) — \r) is a resonant 
dipole-dipole interaction scaling as n 4 . The ratio is thus 
maximized by keeping n not too large. A possible choice 
with Cs atoms is \s) = |60p 3 / 2 ) and |r) = |60si/ 2 ). Us- 
ing a square lattice with the ancilla atom at the center 
we will have k — 1 atoms available for register encoding. 
Averaging the \s) — \r) interaction over the lattice [24] 
we find the error for each full Grover step (controlled 
sign plus inversion about the mean) is (0.04,0.14,0.20) 
for k — 1 = (8,15,24). Although these intrinsic errors 
are a few times larger than for the individual addressing 
protocol, this may be partially compensated by the tech- 
nical advantage that global addressing requires 2/k fewer 
Rydberg pulses. 

The error scaling for the sub-register architecture pre- 



sented in Sec. |IIIB| can be estimated as follows. Denoting 
the error in the simultaneous excitation approach with k 
bits as E(k) the sub-register architecture error for each 



Grover iteration is E SY = n s E(k s ) + (n s — l)E a where 
k = n s k s and E a is the additional error associated with 
each of the ancilla comparison operations. These ancilla 
comparison operations rely on interaction gates, not the 
blockade gates we use elsewhere. Inspection of Fig. 15 
in [5] shows that using, e.g. Rubidium 150s states will 
allow these errors to be less than 0.01 at distances out 
to 100 //m. The physical separation of two centrally po- 
sitioned ancilla atoms in neighboring sub-registers will 
be roughly dy/k~ s so for d < 10 fim and k s < 25 we are 
working with a distance that is under 50 jam. Refer- 
ring to Table [j] we see that the error overhead associated 
with the ancilla comparison steps is therefore negligible 
for k > 8. Since the simultaneous addressing error grows 
faster than linearly with k the sub-register architecture 
can be advantageous and leads to the potential for im- 
plementation on rather large search problems. 

Although the errors shown in the Table are too large 
for standard circuit model fault tolerant operation they 
may be adequate for Grover search, where the idea is 
to increase the amplitude of a desired component of the 
wavefunction. An error in one step may not lead to the 
desired amplification, or may even push the state vector 
in the wrong direction. Nevertheless, as long as the aver- 
age error is sufficiently small we expect the search proce- 
dure to proceed towards the sought after solution. It has 
been shown in [26] that a quadratic speedup is still pos- 
sible provided the oracle phase error per step is bounded 
by TV -1 / 4 . For lar ger errors the quantum speedup will be 
less than quadratic. Comparison of the listed errors with 
the last column in Table [I] shows that a full quadratic 
speedup is possible for problems as large as N = 65536 
with sequential addressing. 

We note that the sequential addressing scheme gener- 
ally has lower errors than simultaneous addressing, with 
or without sub-registers. This conclusion is based on ig- 
noring technical errors associated with imperfect laser 
pulses. Since the simultaneous schemes require a factor 
of k/2 fewer laser pulses they may nevertheless be ad- 
vantageous. A reliable comparison of the approaches, as 
well as an accurate prediction of the possible speedup, 
will depend on implementation details and is beyond the 
scope of this work. 



V. CONCLUSION 

In this paper we have suggested to use the particular 
multi-atom possibilities of the Rydberg blockade inter- 
action to implement the Grover search algorithm using 
multi-bit quantum gate operations rather than the con- 
ventional one- and two-qubit circuit model. The Rydberg 
blockade mechanism "takes the atoms to a third level" 
(the Rydberg state) and hence leaves the qubit and quan- 
tum circuit paradigm for an extended number of opera- 
tions during gate operation. The construction of multi- 
bit operations and entire algorithms without the use of 
predefined universal gates is not viable or meaningful as 
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N 


k-1 


k 


sequential 
addressing 


simultaneous 
addressing 


simultaneous with 
sub-registers 


quadratic speedup 
limit N~ 1/4 


256 


8 






.08 




.25 


512 




9 


.004 






.21 


32768 


15 






.20 




.074 


65536 




16 


.015 




.16 (k s = 8,n s = 2) 


.063 


16777216 


24 






.28 


.24 (k s = 8,n s = 3) 


.016 



TABLE I: Errors per Grover iteration step for the different architectural approaches described in the text. 



a general approach to quantum computing as it merely 
represents a transfer of the computational complexity, 
which should ideally be taken care of by the quantum 
processor, to the theoretical design of the gates. In the 
present case, however, the experimental strategy is so 
straightforward that simple schemes are readily derived. 

We focussed on the Grover algorithm, showing first 
that it can be decomposed into sign changes conditioned 
on all register qubit states, and secondly that these op- 
erations can be accomplished by tt pulse excitations of 
the individual atoms. A general fc-qubit unitary opera- 
tion is described by a 2 k x 2 k matrix, and the theoretical 
lower bound on the number of C — NOT gates needed, 
together with arbitrary single-qubit gates, to form an 
arbitary unitary operation scales as A k [27] . A practi- 
cal, constructive protocol reaching this limit is presented, 
e.g., in [28] . In comparison, a quantum circuit design 
has been made which decomposes the Grover step into 
49/c — 149 elementary one- and two-bit gates (k > 3) [29] . 
Our algorithm requires \fN repetitions of the Grover step 



(controlled sign and an inversion about-the-mean) , and 
we showed that in physical situations where the atoms 
must be excited sequentially this step requires 4fc single 
atom 7r pulses, while for level schemes offering both inter- 
acting and non-interacting atoms in the Rydberg states, 
they can be accomplished in the duration of only 8 sin- 
gle atom 7r-pulses. Furthermore, the estimated fidelity 
benefits from the multi-bit gate operations and is sug- 
gestive that the Grover algorithm can be successfully im- 
plemented on a quite large register. 

We are convinced that our ideas for multi-bit gates, 
which go just a little bit beyond the qubit and circuit- 
model paradigm may be pursued much further, and al- 
though they may not change the speed-up offered by 
quantum computers in a conceptual way, they may of- 
fer practical and also quite efficient ways to carry out 
quantum computing protocols in practice. 

This research was supported by the EU integrated 
project AQUTE, the IARPA MQCO program, and 
DARPA and NSF award PHY-0969883. 
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